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ABSTRACT 

Particle-in-cell/Monte  Carlo  (PIC/MC)  simulations  of  capacitively  coupled  radio-frequency  (RF)  glow 
discharges  were  carried  out  for  low  pressure  CH4/FI2  plasmas.  The  present  computational  scheme  considers  the 
motions  and  collisions  of  both  neutral  and  charged  particles.  The  CH4/H2  plasma  is  modeled  by  combining  a  one¬ 
dimensional  PIC/MC  method  with  a  polyatomic  particle  collision  scheme.  The  model  considers  the  motions  of 

+ 

neutrals  CFI^,  CFIj ,  C^FIj,  and  FI,  positive  ions  (CFl^  )  and  electrons.  Space  and  time  dependent  results  show  that 
ionization  rate  is  high  at  the  sheath  region.  Deposition  rate  of  carbon  film  is  calculated  by  sampling  impinging 
particles  at  the  powered  electrode.  When  CH4  flow  rate  is  decreased  and  the  H2  flow  rate  is  increased,  the  gas  density 
in  the  chamber  and  deposition  rate  decreased  significantly  for  carbon  containing  radicals.  On  the  other  hand,  when 
CH4  flow  rate  is  fixed  and  H  flow  rate  is  increased,  deposition  rate  again  decreased,  but  the  gas  density  in  the 
chamber  did  not  change  much. 

INTRODUCTION 

Diamond  like  carbon  (DLC)  films  deposited  by  various  glow-discharge  processes  have  properties  much  like 
those  of  diamond.  DLC  films  are  extremely  hard,  chemically  inert,  electrically  insulating,  and  good  thermal 
conductors*.  Many  application  areas  have  been  suggested  such  as  spacecrafts  and  heat  sinks  in  electronic 
packaging.  The  role  of  hydrogen  in  vapor-activated  deposition  of  diamond  has  been  the  focus  of  much  attention*'  ^ 
because  only  after  excess  hydrogen  was  used  in  these  processes,  that  mm/h  growth  rates  were  possible.  The  role  of 
hydrogen  is  significant  in  both  gas  phase  chemistry  and  surface  chemistry  for  these  deposition  processes.  In  gas 
phase  chemistry,  atomic  hydrogen  promotes  the  gas  phase  production  of  growth  precursors,  and  suppresses  the 
formation  of  polyatomic  hydrocarbon  species  (PAHs)  by  molecular  hydrogen^.  PAHs  may  be  precursors  to 
unwanted  sp^  carbon  during  the  vapor-deposited  DLC  growth.  Regarding  surface  chemistry,  several  factors  such  as 
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etching,  covering  of  the  dangling  bonds  of  surface  carbons,  removing  hydrogen  from  surface  C-H  bonds  (to  create 
active  growth  sites)  are  affected  by  atomic  hydrogen. 

The  role  of  hydrogen  in  DLC  deposition  process  in  pure  methane  plasma  is  not  clearly  understood'.  Most 
of  the  reported  models  for  DLC  deposition  process  are  for  pure  CH4  plasma ’  These  models  focused  on  gas  phase 
chemistry  and  surface  chemistry  rather  than  discharge  kinetics.  Rate  equations  with  gas  phase  chemistry  or  plug  flow 
model  (simplified  fluid  model)  were  used  to  predict  the  number  densities,  etc.  Electron  energy,  which  determines 
ionization  and  dissociation  rates,  was  obtained  either  by  experimental  data  or  via  Monte  Carlo  simulation  of  electron 
kinetics.  Depositions  were  predicted  using  sticking  coefficients.  A  self-consistent  fluid  model  for  CH4  plasma,  using 
Boltzmann  equation  solver  for  DC  field  calculation  of  the  electron  energy  distribution  function  was  later  developed^. 
The  radical  densities  for  CH3  and  CH2  were  calculated,  however,  these  profiles  did  not  compare  well  with  experimental 
data.  Hybrid  models  have  also  been  developed  and  applied  to  RF  plasma  with  gas  mixtures"*’ ". 

In  this  study,  we  applied  the  PIC/MC  method  to  simulate  CH4/H2  plasmas  in  a  capacitively  coupled  RF 

reactor.  The  model  considers  the  motions  of  CH4,  CH4''“,  CH3  ,  C2H5,  H2,  H  and  electrons.  Both  flow  rates  and  gas 
compositions  or  pressure  were  varied.  The  results  are  analyzed  by  comparing  the  power  consumed,  current  flow,  and 
the  electron  energy  and  density  (at  the  center  of  the  plasma)  as  the  gas  compositions  are  varied.  Ionization  and 
dissociation  rates  and  the  species  densities  are  computed  as  functions  of  the  gas  composition  and  flow  rates.  The 
energy  of  the  impinging  ions  to  the  cathode,  and  deposition  rates  are  also  obtained.  The  developed  computational 
model  will  allow  us  to  predict  the  carbon  film  deposition  rate  and  its  quality  when  different  CH4/H2  compositions  are 
fed  into  low  pressure  plasma  assisted  CVD  reactors. 

MODEL  DESCRIPTION 

The  models  used  here  are  an  extension  of  the  one  dimensional  self-consis  tent  particle  model  for  capacitively 
coupled  RF  glow  discharge  of  pure  A  parallel  plate  configuration  of  the  reactor  is  considered  where  the  two 

electrodes  have  large  cross-sectional  areas  compared  to  their  separation  distance,  so  that  a  one-dimensional 
approximation  is  valid  for  characterizing  the  discharge.  The  model  considers  motion  and  collisions  of  neutrals  as  well 
as  of  charged  particles.  To  simulate  the  CH4  /H2  plasma,  we  combine  the  particle  simulation  model  with  neutral  motion 
and  a  polyatomic  particle  collision  model  which  was  developed  in  the  rarified  gas  dynamics  area'^’  For  stable 
neutral/radical  and  radical/radical  collisions,  energy  after  collision  is  determined  by  the  Borgnakke  and  Larsen"* 
model.  In  this  model,  both  the  relative  transnational  and  the  internal  energies  are  assumed  to  be  distributed 
according  to  the  respective  equilibrium  distributions  and,  at  each  collision,  new  values  are  sampled  at  random  from 
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these  distributions,  subject  to  the  condition  that  the  total  energy  is  conserved.  Hence  vibrational-vibrational  and 
vibrational-translational  energy  exchanges  are  considered  for  the  neutral/neutral  collisions. 

The  model  uses  polyatomic  particle  collision  scheme,  which  has  been  developed  in  the  non-plasma  area. 
Gas  phase  chemistry  is  considered  by  simplified  reactions*^  to  reduce  the  number  of  species  in  the  particle  model.  In 
our  particle  simulation,  we  simplified  the  chemistry  model  of  Kline  et  al.''  Electron-electron,  electron-ion,  and  ion-ion 
collisions  are  neglected  in  the  present  study,  as  the  probability  of  such  collisions  are  small  in  weakly  ionized  plasma. 
For  neutral-electron  collisions,  elastic  scattering,  excitation,  and  ionization  collisions  are  considered.  In  each  collision 
process  total  momentum  is  conserved.  Total  energy  of  collision  pairs  in  elastic  collision  is  preserved,  while  energy  is 
lost  in  excitation  and  ionization  processes.  For  the  discharge  kinetics,  elastic  scattering  between  neutrals  are 
considered  and  electron  impact  to  CH4  and  Fk  are  considered.  Hydrogen  ion  is  neglected  because  the  density  is 
negligible  compared  to  that  of  CH^'*'  in  the  simulated  range  of  parameters.  Following  gas-phase  reactions  are 
considered  for  the  plasma: 

CH4 -he"  ^CH3 -hH+e"  CH3+HH.CH4 

CH4 -he' ^CH^ -he' +e'  CH3 -hCH3  ^C2H5 +H 

H2+e"^H  +  H  +  e‘  H+H^H2 

Surface  reactions  considered  in  the  present  analysis  are  as  follows: 

CH 4  — >  C  (deposit)  +  2H 2  (gas)  CH  3  — >  CH  (deposit)  +  H2  (gas) 

C2H3  — ^C2H (deposit)  +  2H 2  (gas)  H  — >  absorbed  or  H 2  (gas) 

Further  details  of  the  gas-phase  and  surface  reactions  considered  can  be  found  in  Nagayama  et  al.*^.  The  depositing 
species  are  CH/  and  CH3  and  QHs  radicals.  The  deposition  of  hydrogen  radical  is  neglected  in  this  study.  Film 
growth  rate  is  predicted  by  sampling  depositing  particles  on  the  surface  using  sticking  probability.  Sticking 
probabilities  are  assumed  to  be  1.0  for  the  ions  and  0.1  for  the  radicals.  Although  the  model  does  not  consider  a 
detailed  surface  model,  it  provides  useful  information  for  plasma  physics  and  can  identify  major  species  impinging  to 
the  substrate. 

For  the  present  PIC/MC  scheme,  the  reaction  rates  depend  on  the  respective  cross  sections.  In  general,  rate 
constant  for  any  reaction  is  related  to  the  cross  section  for  that  reaction  as 

°°  1 2E 

K=  Jf(E)J - cy(E)dE 

0  V  m 
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where  E  is  the  total  energy  of  the  collision  pair,  f(E)  is  the  energy  distribution  function,  m  is  the  reduced  mass'*  of  the 
collision  pair  and  C5(E)  is  the  cross  section.  For  electron  impact  ionization  or  dissociation,  cross  section  area  is 
expressed  as  a  function  of  electron  energy'*,  which  is  directly  used  in  PIC/MC  simulations.  No  electron-impact 
excitation  reactions  were  considered  in  this  study.  For  neutral  stable/radical  and  radical/radical  reactions,  rate 
constant  must  be  transformed  to  cross  section  o  as  PIC/MC  method  requires  cross  section  area  as  input  data. 
Assuming  hard  sphere  model'*,  the  above  equation  can  be  solved  for  ct(E)  as 

In  order  to  get  stable  results  with  a  relatively  small  number  of  particles  (as  used  in  this  study),  a  time¬ 
averaging  technique  was  used  to  determine  the  electric  field.  Based  on  the  fact  that  RF  plasma  is  a  periodic 
phenomena,  new  electric  field  for  a  given  time  phase  at  the  end  of  the  n*  cycle  is  calculated  using  averaged  density 
data  of  the  previous  n/2  cycles  at  the  same  time  phase.  In  addition  weighting  factors  were  used  to  represent  the 
charged  particles  and  the  neutrals.  Further  details  can  be  found  in  Nagayama  et  al.'^. 

Calculation  conditions  for  the  CHi/ft  plasma 

The  discharge  is  powered  by  RF  source  at  13.56  MHz.  Eighty  uniformly  distributed  grid  points  were  used 
for  the  distance  (2.5  cm)  between  the  electrodes.  1,200  particles  of  each  type  are  used.  The  gas  is  supplied  from  one 
electrode  (anode).  CH^  and  the  neutral  radicals  produced  by  the  reaction  are  partially  eliminated  under  the  condition 
that  total  pressure  inside  the  chamber  is  held  constant.  The  secondary  electron  emission  coefficient  was  set  to  zero 
for  all  calculations  reported  here.  In  surface  reactions,  C-H  or  C2-H  is  absorbed  to  form  film  and  hydrogen  is  released. 
Accordingly,  the  depositing  particles  are  eliminated  from  the  calculation  and  hydrogen  particles  are  added. 
Anisotropic  angular  scattering  is  incorporated  with  the  relative  differential  cross  section  for  all  electron  and  neutral 
collisions.  Electron  wall  reflection  ratio  is  set  to  0.25.  The  calculations  are  carried  out  for  1200  RF  cycles  and  the  data 
are  averaged  between  600th  -  1200th  cycles  for  calculating  the  plasma  properties.  Time  steps  considered  were 
between  0.001  and  1.0  ns  for  a  charged  particle  depending  on  its  velocity.  Larger  time  steps  were  considered  for  the 
neutrals 

For  all  calculations,  the  driving  voltage  (Vrf)  is  set  to  200  V,  and  the  DC  bias  voltage  (Vdc  )  is  equal  to  zero. 
Flow  rate  and  gas  compositions  or  pressure  were  changed.  Two  cases  of  simulations  were  carried  out. 

Case  1:  Total  flow  rate  is  fixed  to  24  seem,  and  supply  gas  composition  CH4/(CH4  +  H2)  ratio  is  varied  keeping 

the  total  pressure  at  400mTorr. 
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Case  2:  CH4  flow  rate  is  fixed  to  24  seem,  and  H2flow  rate  is  changed  keeping  the  CH4  partial  pressure  to 
100  mTorr. 

RESULTS  AND  DISCUSSION 

Case  1 :  Fixed  total  flow  rate  with  varying  composition 

Simulations  were  carried  out  at  different  H2  partial  pressures  (300,  200,  and  0  mTorr),  while  total  pressure  is  kept 
constant  (400  mTorr).  In  other  words,  the  gas  composition  is  changed  while  total  pressure  is  kept  constant.  The  simulated 
power  consumption  and  current  in  the  plasma  (as  a  function  of  the  gas  composition)  are  shown  in  Fig.  1-a. 
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Figure  1-a.  Power  and  displacement  current  as  function  of  CH4/(CH4  +  H2) 


Figure  1-h.  Electron  energy  and  electron  density  (at  the  center  of  the  discharge)  as  function  of  CH4/(CH4  +  H2) 


Note  that  small  values  of  CFl4/(CH4  +  H)  value  indicates  that  the  gas  composition  is  highly  diluted  by 
hydrogen.  Power  consumed  in  the  reactor  does  not  depend  on  the  gas  composition  in  the  simulated  range.  Total 
current  tends  to  increases  as  CH4  ratio  increases,  which  may  be  related  to  the  increase  of  plasma  density.  Fig.  1-b 
shows  the  electron  density  and  electron  energy  values  at  the  center  of  the  discharge  as  a  function  of  the  above  ratio. 

Electron  density  increase  as  CH4  ratio  increases  because  the  CH4  ionization  is  the  major  source  of  electrons. 
Electron  energy  is  between  3.4  -  3.8  eV  in  the  simulated  range.  The  electron-impact  ionization  and  dissociation  rates 
are  shown  in  Fig.  1-c  as  function  of  the  inlet  gas  composition  ratio. 
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Figure  1-c.  Computed  (electron-impact)  reaction  rates  as  a  function  of  CH4/(CH4  +  H2) 

CH^  dissociation  rate  is  about  ten  times  higher  than  CH^  ionization  rate  and  the  variations  are  small  in  the 
log  scale.  Although  dissociation  rate  is  negligible  in  pure  CH4  plasma,  H2  dissociation  rises  to  the  same  level  as 
CH4  ionization  rate  at  highly  diluted  condition.  Although  we  neglected  FI  ionization  process  in  this  model,  the 
assumption  is  reasonable  because  such  rate  will  be  lower  than  H2  dissociation  and  will  be  negligible  compared  to 
CH4  ionization  and  dissociation  in  the  simulated  range. 

Fig.  2  shows  the  radical  and  neutral  densities  as  functions  of  the  CFl4/(CFl4  +  FI2)  ratio.  CFI4  increases  and  FI2 
decreases  as  the  gas  composition  is  changed.  H  density  does  not  change  much  although  density  increases 
significantly.  The  reason  may  be  that  FI  is  generated  from  both  CH4  and  H2  dissociation.  CH^  density 


Figure  2.  Computed  neutral  density  profiles  as  functions  of  CFl4/(CFl4  -I-  H2) 
increases  as  and  CH4  dissociation  increases.  C^Hj  density  also  increases  as  CH^  density  increases. 

The  energy  of  ions  impinging  to  the  cathode  tends  to  decrease  as  CFl4  density  increases  and  charge 
exchange  collisions  increase.  It  is  noted  that  the  present  model  does  not  consider  the  interaction  between  CH4  ion 
and  H2.  Fig.  3  shows  the  deposition  rates  for  the  three  depositing  species  as  function  of  the  CFl4/(CH4  +  H2)  ratio.  CH^ 
and  deposition  rates  increase  as  CH4  density  increases  C^Hj  represents  PAHs  (polyatomic  hydrocarbon)  and 
high  may  cause  polymer  like  films.  On  the  other  hand  deposition  rate  is  constant  because  ionization  rate 
does  not  change. 
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Figure  3.  Deposition  rates  as  function  of  functions  of  CH4/(CH4  +  H2) 

Case  2:  CH^  flow  rate  is  fixed  while  hvdroeen  is  added 

Simulations  were  also  carried  out  at  different  Ffc  flow  rates  (0,  48,  96  seem),  while  CH4  flow  rate  is  kept 
constant  at  24  seem.  The  CH4  partial  pressure  is  thus  indirectly  fixed  at  100  mTorr.  The  power  consumption  and 
current  are  high  for  high  H  addition.  It  means  power  is  mostly  consumed  by  electron-hydrogen  collision. 
Although  FI2  density  increases  as  the  dilution  becomes  high,  the  density  levels  of  most  species  do  not  change  as  in 
the  previous  case.  The  ion  energy  is  between  24  -  28  eV  and  does  not  show  any  variation  with  hydrogen  addition.  As 

noted  before,  no  interactions  between  and  H^,  were  considered  here. 

Fig.  4  shows  the  spatial  profiles  of  time  averaged  rate  coefficients  for  different  values  of  pressure  as  CH4 
pressure  is  held  constant  at  100  mTorr. 


Figure  4.  Time  averaged  rate  coefficient  profiles  for  electron-impact  reactions  with  varying  H2  addition 

With  pure  CFI4,  dissociation  rate  is  large  and  dissociation  is  negligible  (Fig.  4-a).  The  CFL*  dissociation  rate  is 
high  at  the  center  of  the  plasma  that  is  typical  at  low  pressure  (100  mTorr).  As  the  pressure  is  increased  (Fig.  4-b,  c), 
dissociation  is  increased.  It  is  interesting  to  note  that  by  adding  H,  the  spatial  profile  of  the  rate  coefficients  become  bow 
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shaped,  which  is  similar  to  the  case  of  pure  gas  discharges  where  pressure  is  increased.  At  pressure  of  400  mTorr  (80% 
of  the  total  flow),  dissociation  rate  is  comparable  to  CH,  ionization  rate  but  still  about  ten  times  smaller  than  CH4 
dissociation  rate.  dissociation  and  CH4  ionization  rates  are  high  in  the  sheath  region  while  CH4  dissociation  rate  is 

uniform  in  the  plasma  bulk  at  this  pressure  (Fig.  4-c).  The  reason  may  be  due  to  the  fact  that  in  the  former  two  processes, 
high  energy  (threshold  energy  is  13  eV  for  dissociation  and  14  eV  for  CH4  ionization)  is  needed  compared  to  CH4 

dissociation  (threshold  energy  is  10  eV). 

CONCLUSIONS 

Particle-in-cell/Monte  Carlo  simulations  of  capacitively  coupled  RF  (radio-frequency)  glow  discharges  were 
carried  out  for  low  pressure  CH^  /H^  plasmas.  The  model  employs  a  one-dimensional  scheme  and  considers  polyatomic 

particle  motion  and  collision.  Calculations  were  carried  out  for  two  cases;  (a)  total  flow  rate  is  fixed  while  gas  composition  or 
CFl4/(CH4  -I-  H2)  ratio  is  varied,  and  (b)  CH4  flow  rate  is  fixed  while  H2flow  rate  is  changed.  When  CH4  flow  rate  is  decreased, 

4- 

the  gas  density  in  the  chamber  and  deposition  rate  decrease  significantly  for  carbon  containing  radicals  while  CFl^  density 
does  not  change  much.  On  the  other  hand,  when  CH4  flow  rate  is  fixed  and  H2  flow  rate  is  increased,  and  deposition  rate 
decrease  similarly,  but  the  gas  density  in  the  chamber  does  not  change  much  (for  the  range  of  parameters  considered).  For 
both  cases,  FI2  dissociation  rate  is  comparable  to  CH4  ionization  rate  when  H2  partial  pressure  is  about  75  ~  80%  of  the  total 
gas  pressure.  Flowever,  dissociation  rate  is  still  very  small  compared  to  CH4  dissociation  rate.  When  H  flow  rate 
increases  with  fixed  CFI4  flow  rate,  lower  CH^  and  C^Hj  deposition  rates  occur,  but  the  overall  effect  is  not  large  (H2  flow 

rates  between  0  to  80%  of  the  total  flow).  The  major  reason  is  that  electron  impact  dissociation  cross  section  of  hydrogen 
is  smaller  than  that  of  CH4.  With  the  addition  of  H2,  carbon-containing  gas  decreases  significantly  and  deposition  rate  is 
decreased  especially  for  C^H..  while  that  for  does  not  change  much.  This  may  improve  the  film  quality  suppressing 
the  PAHs  formation.  Adding  also  changes  the  rate  coefficient  spatial  profiles  to  bow  shape,  which  is  similar  to  the  case 
of  pure  gas  discharges  when  pressure  is  increased. 


REFERENCES 

1 .  K.  E.  Spear  and  J.  P.  Dismukes,  Synthetic  Diamond:  Emerging  CVD  Science  and  technology.  John  Wiley  & 
Sons,  Inc.  (1994) 

2.  M.  Frenklach  ,  J.  Appl.  Phys.  65  (12),  5142  (1989) 

3.  Y.  Catherine  and  P.  Coudrec,  Thin  Solid  Films  144, 165  (1986) 

4.  L.  E.  Kline,  W.  D.  Partlow,  and  W.  E.  Biee,  J.  Appl.  Phys.  65  (1),  70  (1989) 

5.  A.  Rhallabi  and  Y.  Catherine,  IEEE  Trans.  Plasma  Sci.  19  (2),  270  (1 991 ) 

6.  E.  Gogolides,  D.  Mary,  A.  Rhallabi,  G.  Turban,  Jpn.  J.  Appl.  Phys.  34  261  (1995) 

7.  A.  von  Keudell  and  W.  Moller,  J.  Appl.  Phys.  75  (12),  7718  (1994) 

8.  K.  Tachibana ,  M.  Nishida ,  H.  Haiima ,  Y.  Urano,  J.  Phys.  D:  Appl. .  Phys.,  1 7, 1727  (1984) 

9.  E.  Gogolides,  D.  Mary,  A.  Rhallabi,  G.  Turban,  Jpn.  J.  Appl.  Phys.  34  261  (1995) 

10.  T.  J.  Sommerer  and  M.  J.  Kushner,  J.  Appl.  Phys.  71  1654  (1992) 

11.  N.  Sato  and  H.  Tagashira,  IEEE  Tran.  Plasma  Sci.  19  102  (1991) 

12.  K.  Nagayama,  B.  Farouk  and  Y.H.  Lee,  IEEE  Tran.  Plasma  Sci.  26  (2),  125  (1998) 

13.  G.  A.  Bird ,  Molecular  Gas  Dynamics  and  The  Direct  Simulation  of  Gas  Flow,  Oxford  University  Press, 

New  York,  NY,  1994 

14.  C.  Borgnakke  and  P.  S.  Larsen,  J.  Comp.  Phys.  18, 405,  (1975) 

1 5.  D.  J.  Economou,  T.  J.  Bartel,  R.  S.  Wise  and  D.  P.  Lymberopoulos,  IEEE  Trans.  Plasma  Sci.  23, 581  (1995) 

16.  W.  L.  Morgan,  Plasma  Chem.  Plasma  Process.,  12  (4),  477,  (1992) 


237 


